Molecular spintronics using noncollinear magnetic molecules 
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We investigate the spin transport through strongly anisotropic noncollinear magnetic molecules 
and find that the noncollinear magnetization acts as a spin-switching device for the current. More- 
over, spin currents are shown to offer a viable route to selectively prepare the molecular device 
in one of two degenerate noncollinear magnetic states. Spin-currents can be also used to create a 
non-zero density of toroidal magnetization in a recently characterized Dy3 noncollinear magnet. 



One of the most ambitious directions in the quest for 
the ultimate miniaturization of electronic devices is rep- 
resented by molecular spintronics [TJ |2] . Molecular nano- 
magnets are particularly promising for nanospintronics, 
especially in relation to the quest for magnetic molecu- 
lar qubits [3 , since transport experiments [4j [5j [6] have 
shown a strong interplay between the current and the 
magnetic states of the molecules. To date, all theoretical 
investigations on molecular spintronics have addressed 
systems whose magnetism is only weakly anisotropic, 
thus exploring systems whose magnetization aligns along 
a single anisotropy axis (collinear magnetism) [TJ [H |9j 
[TO] . The noncollinear regime of molecular magnetism, 
arising when the on-site magnetic anisotropy of single 
metal ions is one of the dominant energy scales, has only 
been explored very recently [11 j. In the noncollinear 
regime magnetic molecules can be prepared in degenerate 
states characterized by non-dipolar magnetic moments, 
such as the recently predicted [11] and found [12j [13] 
toroidal (or anapole) moment in molecular wheels. 

There are two main arguments that make spintron- 
ics of noncollinear magnets of great interest. The first 
follows from studies on spin-transport through meso- 
scopic rings with noncollinear internal magnetic fields, 
which have been predicted to produce spin-switching ef- 
fects [14 J. The size of noncollinear molecular scatterers 
is expected to be more favorable to overcome dephas- 
ing, and lead to the observation of coherence and spin- 
switching effects. The second argument is related to the 
use of noncollinear states to implement molecular qubits. 
On the one hand, molecular spin-qubits can easily be ad- 
dressed via a magnetic field [15] , although intermolecular 
dipolar interactions lead to short dephasing times [T6] . 
On the other hand, intermolecular interactions between 
non-dipolar states are weak [TT] . thus decoherence times 
longer, although these states cannot be addressed via uni- 
form fields. Spintronics might offer a promising strategy 
to address noncollinear protected molecular qubits. 

In this Letter we investigate spin transport through 
molecular noncollinear magnetic states, and provide ev- 
idence that these systems do offer strategies to (i) im- 
plement quantum-interference molecular devices capable 
of reversing the polarization of an injected spin-current, 
and to (ii) selectively populate specific noncollinear mag- 




FIG. 1: (color online) The two molecular spintronics setups 
considered in this work. Top, a noncollinear antiferromag- 
netic dimer with on-site spin s — 3/2 and coplanar on-site 
ZFS easy-axes, tilted with respect to the perpendicular to 
the inter- metal distance by ±30°. Bottom, a 3-center anti- 
ferromagnetic molecular wheel with on-site spin s — 3/2 and 
coplanar on-site ZFS easy- axes, arranged tangentially to the 
wheel's circumference (cf Ref. [Ill 112] ). 

netic states. The most relevant transport regime has 
been shown to be the Coulomb blockade (CB) [6 . The 
lowest lying states of a nanomagnet with n unpaired elec- 
trons well localized on N metal centers with local spin s 
is well described by the Hamiltonian: 

H n = -J^Si-Sj+Dj^sh (1) 

(ij) i 

consisting of the isotropic Heisenberg exchange coupling 
between nearest neighbors with strength J, and easy- 
axis zero-field splitting (ZFS) on-site (strength D with 
D < 0). Note that the spin operator Si has z-component 
Si yZ parallel to the local ZFS axis. Whereas previous 
investigations [10] only considered the collinear weak- 
anisotropy regime, here we introduce two key-ingredients 
for noncollinearity [TT]: (i) \D\ » \J\ (ii) on-site easy- 
axes not parallel to each other. In this work we will 
explore spin-transport for a dimer (Fig. 1, top) and a 
three-centers molecular wheel (Fig. 1, bottom). 

When connected to source, drain and gate electrodes, 
under bias voltage Vb and gate voltage Vq, the molecule 
will become charged. The migrating electron will be con- 
secutively accommodated at different metal sites [17], de- 
scribed here by a set of N atomic orbitals localized on the 
metal centers. The molecular Hamiltonian for a charged 
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state with an excess of Q electrons with respect to the 
isolated nanomagnet is given by: 

n n 

H n+Q = H n + (e - eQV G ) ^2 ^2 n P a 

n n 

(pq) cr p 

n n n 

+ J ^^^J2 g p' a p^ c U c Pi3^ ( 2 ) 

p a (3 

where e is the energy of the localized orbit als, c^ po . are 
creation operators for the on-site spin-orbit als, n pa = 
c \)cr c pcn ^ is a hopping parameter between centers, U is 
the Coulomb repulsion between two electrons on the same 
center, a p are Pauli matrices associated to an electronic 
spin injected on site p, and Jh is the Hund's rule cou- 
pling between the spin of the excess electron on site p, and 
the spin moment s p on that center (Jh < 0). Here we 
confine ourselves with the region around the first CB di- 
amond, where only singly-charged states are relevant to- 
gether with the neutral ones. This is formally achieved by 
setting U — > oo. Finally, a tunneling Hamiltonian H m i x 
between electrodes and device is introduced in the usual 
manner pj [H [10] , with tunneling amplitudes estimated 
to be at most 0.3 cm -1 [6j[9]. Given the weak molecule- 
lead coupling, the transition rates W between molecule 
and contacts are calculated with the Fermi golden rule us- 
ing H mix , assuming a Fermi-Dirac distribution in the two 
leads, kept at different chemical potential Hl—^r = eVs- 
Next, using the rates W, a master equation for the non- 
equilibrium populations of charged and neutral states of 
the device is set up and solved in the steady-state regime. 
The resulting populations are used to compute the input 
(l\ — I L ) and the output (J^ — ijfc) spin currents. 

Let us first consider the dimer molecule with local spin 
s = 3/2 (e.g., a Co(II) dimer). We assume co-planar local 
ZFS axes, forming angles — ±30° with the perpendic- 
ular (z-direction) to the Co-Co bond, (Fig. 1, top). We 
choose here D — 5 J, with antiferromagnetic isotropic 
exchange coupling J = — 50 cm -1 . The dominant en- 
ergy scale in ([2| is the on-site Hund coupling exchange 
Jh, chosen here as Jh = 4D « 0.1eV\ Moreover, the 
source contact is assumed to be ferromagnetic, the drain 
non-magnetic. The spin-polarization axis for the ferro- 
magnetic source is coplanar to the ZFS axes, and parallel 
to z (Fig. 1). The Heisenberg states of lower energy of 
the Co-dimer can be described in terms of almost pure 
noncollinear Ising states | m\ 7712), where mi is the pro- 
jection of the local spin moment s along the tilted easy 
axis [11 . This is verified by decomposing the Heisen- 
berg wavefunction into the noncollinear Ising basis. For 
the present choice of parameters, denoting "+" and "— " 
the on-site spin-states |±3/2), the ground state is quasi- 
degenerate and corresponds to the in-phase and out-of- 
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FIG. 2: (color online) Input (blue bullet datapoints) and out- 
put (purple square datapoints) spin-current vs. voltage curves 
(left column), and, (right column) non-equilibrium popula- 
tions of neutral (blue bullet for in-phase and purple square 
for out-of-phase superposition of \-\ — ) and | — h)) and charged 
states (green upward-pointing triangles for ^q^ 1 , yellow di- 
amonds for ^q^ 1 ) f° r t ne noncollinear magnetic dimer, with 
the following type of exchange interaction and values of the 
transfer parameters (D — —250 cm -1 ): a) and b) Heisen- 
berg, t — 0.5£>, c) and d) Heisenberg, t = 4.0D, e) and f) 
noncollinear Ising, t — 1.0D. 

phase superposition of the noncollinear Neel states |H — ), 
and I — h). The tunneling gap is about A w 0.38 cm -1 . 

Next, we find the eigenstates of We explore here 
two limiting situations: a weak-transfer limit with t = 
0.5D, and a strong-transfer limit with t = AD = Jh- In 
the weak-transfer limit, we expect the additional electron 
spin to follow "adiabatically" the noncollinear magnetic 
texture of the molecular device, so that the ground state 
of the charged system will be doubly degenerate, and 
dominated by either the |H — ) or the | — h) component 
of the Neel doublet, carrying an additional electron on 
either center, with the spin parallel to the local magnetic 
moment via Hund-rule coupling. These expectations are 
confirmed by full diagonalization, leading to the following 
ground state for the n + 1-magnet: 

^o,! 1 « C T , |+-> IT 0) + C ,x |+-> |0 i) 

l-o.t 1 « C uo |-+) || 0) + Co, T |-+) |0 T) , (3) 
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where the kets || 0) denote the determinant of two spin- 
orbit als centered on the two metals, the first with spin 
|, the second empty. We find |C^o| 2 = |C|,o| 2 « 0.45, 
and |Co,J 2 = | Co,t 1 2 ~ 0-42. On the other hand, in 
the strong-transfer limit, inter-center hopping processes 
have the same rate as Hund-coupling spin-polarization 
processes, so that transport is not expected to be adia- 
batic. This is confirmed by full diagonalization, where 
the ground state is dominated by noncollinear Ising 
states favoring spin-preserving hopping processes, such 
as |++) |T 0) and |++) |0 |). Thus, in this regime the 
overlap with the neutral noncollinear Neel states is very 
small. The voltage Vq is taken large enough to bring in 
resonance the ground states of the neutral and charged 
systems, separated by about leV. Without loss of gen- 
erality, we set the equilibrium chemical potential lying in 
between the ground and first excited state of the neutral 
molecule, and the temperature to T = 0.1K. 

In Fig. 2a (weak-transfer) and 2c (strong-transfer) we 
report the spin current- voltage diagrams obtained for the 
two limits of the hopping parameter. Since the source is 
fully spin-polarized, the input spin current (blue bullets) 
always corresponds to the total charge-current. Interest- 
ingly, in the weak-transfer limit the output spin-current 
(purple square datapoints) has a negative sign: the spin 
polarization of the input current is reversed in the out- 
put non-magnetic electrode, by the noncollinear magnetic 
texture. On the other hand, in the strong-transfer regime 
this spin- switching effect is not observed. These results 
are easily interpreted analyzing the ground state wave- 
functions for the neutral and charged states. In the weak- 
transfer limit, the charged ground state Q is a coher- 
ent state describing the adiabatic hopping of an injected 
electron-spin between the two metals, in which process 
the additional electron always aligns its spin parallel to 
the magnetic polarization of the local metal ion. Thus 
we define this limit as the adiabatic-transport limit, in 
analogy with the findings reported in ref. [14 Although 
the CB-regime is non-coherent, the transition rates en- 
tering the master equation are determined by the over- 
lap amplitudes between the tunneling combinations of 
the ground noncollinear Neel doublet, and the charged 
ground-state doublet (J3|. Due to the full |-spin polar- 
ization of the source, the injected electron on the first 
metal center creates an excess of non-equilibrium popu- 
lation in the state ^q^ 1 (see Fig. 2b, green triangular 
datapoints), which can host an electron with spin- up on 
the first metal. The electron is then coherently trans- 
ported through ^q^ 1 on the second metal center, where, 
as described by the |0 j) component of ^q^ 1 , its spin- 
polarization is reversed. Output tunneling events from 
the second metal center into the drain will thus occur 
more frequently with opposite spin-polarization. Since 
the additional spin-polarized electron collapses the tun- 
neling wavefunction into one of the two Neel states, at 



non-zero bias voltage we note that the ground state and 
first excited tunneling states of the neutral system be- 
come equally populated (see Fig. 2b and 2d). 

We note that the coupling between \-\ — ) || 0) and 
H — ) |0 I) m (HI? which determines the spin-switching 
transport, is triggered by the Hund-Hamiltonian. Im- 
portantly, if the angle is set to zero, i.e. within the 
collinear regime, the Hund- mechanisms leading to the su- 
perposition ([3| are not active. Hence, noncollinearity is 
found to be a crucial ingredient for the realization of the 
spin- switching effect. In the non-adiabatic regime, the 
spin-switch effect is quenched (Fig. 2c) due to the neg- 
ligible presence of spin-switch coherences in the charged 
ground state. In this limit the large hopping integral 
favors spin-preserving hopping processes. The current 
magnitude is also significantly smaller due to the small 
overlap between charged and neutral states (see Fig. 2c). 

It is interesting to investigate the case of exact degener- 
acy between \-\ — ) and | — h), implied by the noncollinear 
Ising exchange Hamiltonian: 

H n = -JiJ2 §i > z ^> z + D J2 S h- ( 4 ) 

(ij) i 

The ground state of the charged system for t = D is 
still doubly degenerate and given by ([H]), with |C^o| 2 = 

lOu I' = \Cuo I 2 = \CoM 2 = 0.41. The full T-spin po- 
larization of the source electrode will favor population- 
transfer processes mainly between the \-\ — ) Neel state 
and the spin-switch excited state ^q^ 1 , via the Hund- 
coupling mechanism. However, now only half of the in- 
put j-current is converted to j-current, so that the cur- 
rent is non-polarized in the drain (see Fig 2e). This 
can be rationalized in terms of additional population- 
transfer between | — h) and ^q^ 1 ' smce the latter con- 
tains 7% of the Hund-unstable component | — h) |T 0). 
However, the dominant population-transfer process re- 
mains |H — ) — > ^ol 1 ' an d> as seen f rom Fig. 2f, this fact 
has a fundamental consequence: at non-zero bias voltage 
the spin-current causes a net excess of population of one 
of the two degenerate noncollinear Neel states. Thus, the 
neutral system is prepared in the | — h) state. 

Finally, we consider a 3-center molecular wheel with 
local ZFS axes contained in the molecular plane and 
tangential to the wheel's circumference. This system 
is of special interest, being a model for the experimen- 
tally characterized lanthanide wheel Dy3 [18], which has 
been recently shown to have almost tangential on-site 
anisotropy axes, leading to toroidal magnetization [T2] . 
For simplicity, here we consider an analog of this system 
with s = 3/2 on metal sites. The collective states are 
modelled by the noncollinear Ising Hamiltonian Q, with 
ferromagnetic exchange Jj — 25 cm -1 , and easy- axis ZFS 
parameter D = 8 J/. The ground state of the 3- wheel is 
a doubly degenerate Kramer's doublet characterized by 
a toroidal magnetic moment r = ^bR^2 p s Z}P pTJ H2], 
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FIG. 3: (color online) Input (blue bullet datapoints) and out- 
put (purple square datapoints) spin-current vs. voltage curves 
(left column), and, (right column) non-equilibrium popula- 
tions of neutral (blue bullet for | ) and purple square for 

| + + +)) for the 3-center noncollinear magnetic wheel, with 
the following values of the transfer parameters: a) and b) 
t = 0.05L>, c) and d) t = 1.0D (D = -200 cm" 1 ). 

where R is the radius of the triangle, and \±b is the 
Bohr magneton. We denote the two states with | + + +) 

(r = +9/2Rfi B ) and | ) (r = -9/2R{i B ), where the 

first position refers to the atom more strongly bound to 
the ferromagnetic source, and the second position refers 
to the atom bound to the non-magnetic drain. The 
singly-charged system is investigated for J B = 4D, and 
for t = 0.05D (adiabatic transfer) and t = D (strong- 
transfer). The ground state of the singly charged sys- 
tem is always 4- fold degenerate. The present spintronics 
setup (bottom of Fig. 1, spin-polarization axis of the 
source co-planar with the wheel's plane, and perpendic- 
ular to the bond between metal 1 and metal 2) implies 
that only those components of the charged ground state 

overlapping with the | ) toroidal state reported in 

Fig. 1 will be significantly populated, by virtue of the 
Hund's coupling rule. In the adiabatic limit these states 
correspond to spin-switching states, i.e. to states which 
represent coherent hopping from center 1 to center 2, 
with inversion of spin-polarization, and found to be: 

^o,! 1 » I > (ax |t 00) + 6x |0 | 0» 

~ I ) 

x (02 IT 00) + b 2 |0 | 0) + c 2 |00 T) + c 2 |00 I)) 

with |ai| 2 = 0.4, |&i| 2 = 0.38, \a 2 \ 2 = 0.21, \b 2 \ 2 = 0.27, 
and |c 2 | = 0.19. In the non-adiabatic limit the weight 
of the spin- switching components becomes smaller (al- 
though does not vanish), in favor of states representing 



spin-conserving hopping processes. 

In figure 3a (adiabatic) and 3c (nonadiabatic) we re- 
port the spin current-voltage diagrams obtained for the 
two limits of the hopping parameter: as for the dimer 
system, we observe spin switching only in the adia- 
batic limit. However, due to the spin-polarization of 
the source electrode, the population transfer from the 
| ) toroidal neutral state to the charged manifold al- 
ways dominates the non-coherent kinetics, producing an 
excess of population of |+ + +), in both weak and strong 
transfer limits (Fig. 3b and 3d). This demonstrates a 
viable spintronics strategy to prepare a non-zero density 
of toroidal molecular magnetization in the sample. 

In conclusion, we have investigated spin-transport 
through noncollinear magnetic molecules in the sequen- 
tial tunneling (CB) regime. Two fundamental phenom- 
ena are identified here. The first, the spin- switching ef- 
fect, is caused by the action of the noncollinear mag- 
netization on the spin-current. The second, the selec- 
tive population bias of one of the two partners of a non- 
collinear doublet, is determined by the effect of a spin- 
current on the noncollinear states. Non-collinearity is 
found to be the crucial ingredient in these phenomena. 
This work represents the first step into the new domain 
of noncollinear molecular spintronics, expected to have a 
significant impact on the quest for protected molecular 
qubits. 



* Electronic address: lAles san- 

dro.Soncini@chem.kuleuven.be 
[1] A. Rocha et al, Nature Mater. 4, 335 (2005); 
[2] L. Bogani and W. Wernsdorfer, Nature Mater. 7, 179 

(2008). 

[3] G. A. Timco et al, Nature Nanotech. 4, 173 (2009). 

[4] J. Park et al, Nature 417, 722 (2002). 

[5] W. Liang et a/., Nature 417, 725 (2002). 

[6] H. B. Heersche et a/., Phys. Rev. Lett. 96, 206801 (2006); 

M-H. Jo et a/., NanoLett. 6, 2014 (2006). 
[7] F.Elste and C. Timm, Phys. Rev. B 71, 155403 (2005). 
[8] C. Romeike et a/., Phys. Rev. Lett. 96, 196601 (2006). 
[9] G. Gonzales and M.N. Leuenberger, Phys. Rev. Lett. 98, 

256804 (2007). 

[10] M. Misiorny and J. Barnas, Phys. Status Solidi B 246, 
695 (2009). 

[11] A. Soncini and L. F. Chibotaru, Phys. Rev. B 77, 

220406(R) (2008). 
[12] L. F. Chibotaru, L. Ungur and A. Soncini, Angew. Chem. 

Int. Ed. 47, 4126 (2008). 
[13] J. Luzon et a/., Phys. Rev. Lett. 100, 247205 (2008). 
[14] D. Frustaglia, M. Hentschel and K. Richter, Phys. Rev. 

Lett. 87, 256602 (2001). 
[15] F. Troiani et al, Phys. Rev. Lett. 94, 207208 (2005). 
[16] A. Morello et a/., Phys. Rev. Lett. 97, 207206 (2006). 
[17] L.F. Chibotaru et a/., J. Am. Chem. Soc. 125, 12615 

(2003). 

[18] J. Tang et a/., Angew. Chem. Int. Ed. 45, 1729 (2006). 



